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Chapter  1.  SYSTEM  OVERVIEW 


V.  1  Importance  of  Temperature  in  Semiconductor  Measurements 

It  is  well  known  that  most  electrical  properties  of  semiconductor 
materials  are  dependent  in  some  manner  on  temperature.  One  example  is 
the  rate  of  thermal  emission  of  trapped  holes  or  electrons  from 
generation-recombination- trapping  centers.  The  transient  capacitance 
method  can  be  used  to  measure  the  thermal  emission  rate.  This  measure¬ 
ment  is  more  conveniently  accomplished  at  low  temperatures  since  the 
emission  rate  is  smaller  and  thus  the  duration  of  the  capacitance 
transient  is  longer.  [1]  It  is  clear  that  if  accurate  measurements 
of  such  temperature  dependent  properties  are  to  be  obtained,  the 
temperature  must  be  accurately  maintained. 

The  purpose  of  this  thesis  work  is  to  design,  implement,  and  test 
a  computerized  temperature  control  system  which  will  become  part  of  a 
highly  automated  data  acquisition  system.  This  data  acquisition  system 
was  previously  developed  by  the  author  and  other  members  of  the  Solid 
State  Electronics  Laboratory.  Stable  temperatures  in  the  range  from 
85  K  to  300  K  are  obtainable  with  this  system.  In  actual  use,  the  most 
common  input  to  the  controller  will  be  the  step  input.  Therefore,  the 
design  i3  aimed  at  obtaining  a  suitable  step  response. 

1.2  Thesis  Organization 

In  the  remaining  sections  of  this  chapter  a  description  of  the 
physical  apparatus  used  io  construct  the  control  system  is  presented. 

An  approximate  model  cf  the  thermal  portion  of  the  system  is  obtained. 

Two  design  techniques  were  employed  in  the  search  for  a  suitable 
control  a.gorithm  to  be  °xecuted  by  the  computer.  Chapter  2  covers 


the  design  of  a  discrete-time  controller  using  the  z-transform.  In 


Chapter  3  the  system  is  treated  as  an  approximation  to  a  continuous- 
time  system.  The  design  of  a  controller  is  carried  out  using  the 
Laplace  transform.  The  controller  which  gave  the  best  performance  was 
designed  using  Laplace  transform  techniques.  Testing  of  this  version 
of  the  temperature  control  system  is  also  covered  in  Chapter  3.  Chapter 
4  is  the  conclusion  in  which  the  capabilities  and  limitations  of  the 
system  are  summarized. 

1.3  The  Data  Acquisition  System 

The  heart  of  the  automated  data  acquisition  system  mentioned  in 
Section  1.1  is  a  Hewlett-Packard  1000M  minicomputer  system.  Programs 
are  executed  under  control  of  the  Hewlett-Packard  Real  Time  Executive  IV 
(RTE-IV)  operating  system.  In  order  to  gather  data  related  to  the 
properties  of  semiconductor  materials,  various  instruments  and  meters 
must  be  interfaced  to  the  computer.  The  Hewlett-Packard  Multiprogrammer 
Model  694QB  provides  an  easy  and  efficient  solution  to  this  interface 
problem. 

The  purpose  of  the  6940B  Multiprogrammer  is  to  extend  the  capa¬ 
bility  of  one  computer  I/O  channel.  It  is  a  bidirectional  I/O  inter¬ 
face  which  facilitates  communication  between  the  computer  and  several 
input/output  circuit  cards.  Each  6940B  can  contain  up  to  fifteen 
input/output  cards.  In  addition,  up  to  fifteen  6941B  Extender  units, 
each  capable  of  holding  fifteen  cards,  can  be  connected  to  the  6940B. 

In  this  manner  one  computer  I/O  channel  can  be  expanded  to  a  total  of 
240  subchannels.  The  digital  data  sent  to  the  Multiprogramraer  is  used 
to : 

(1)  control  the  mode  of  operation  of  the  system,  for  example 


to  select  whether  data  is  to  be  input  or  output, 

(2)  address  output  cards  and  supply  the  output  quantity  in 


•3 


binary  encoded  form,  or 

(3)  address  input  cards  so  that  binary  encoded  data  may  be 
received  from  the  card. 

Examples  of  the  functions  of  output  cards  are  resistance  programming  of 
power  supplies  and  digital-to-analog  voltage  conversion.  Conversion 
of  an  analog  voltage  to  digital  data  is  an  example  of  the  functions  of 
input  cards.  [2] 

A  block  diagram  of  the  data  acquisition  system  hardware  is  shown 
in  Figure  1.1.  A  detailed  drawing  of  the  probe  which  holds  the  semi¬ 
conductor  device-under-test  and  also  constitutes  the  thermal  portion 
of  the  temperature  control  system  is  shown  in  Figure  1.2.  Power  supplies 
No.  1  and  No.  2  are  programmed  by  the  value  of  resistance  set  on  the 
resistance  programming  cards  in  slots  403  and  404.  The  outputs  of  these 
power  supplies  are  routed  to  the  relay-readback  card  in  slot  406.  Eelay 
closures  on  this  card  select  which  power  supply  voltage  is  applied  to 
the  bias  terminals  of  the  Eoonton  capacitance  meter.  The  polarity  of 
the  bias  voltage  is  also  determined  by  the  relay-readback  card.  The 
power  supply  voltages  are  also  sent  to  one  of  the  voltage  monitor  (A/D) 
cards  in  slots  411  through  413  when  the  appropriate  relay  is  closed. 
Voltage  monitors  001,  010,  and  102  can  measure  voltages  in  the  ranges 
♦1. 0235V  to  -1.0240V  in  0.5mV  steps,  +10. 235V  to  -10.240V  in  5mV  steps, 
and  +102. 35V  to  -102. ^0V  in  50mV  steps  respectively.  Voltage  monitor 
001  may  also  receive  the  analog  output  of  the  Eoonton  capacitance  meter 
through  the  relay  card  in  slot  UQ6.  The  bias  voltage  and  analog  output 
of  the  Eoonton  capacitance  meter  are  multiplexed  to  the  input  of  the 
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Figure  1.1  Elock  diagram  of  the  data  acquisition  system. 
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Data  Precision  3500  digital  voltmeter,  labelled  DVM  2,  by  the  relay- 
readback  card  in  slot  405.  The  maximum  resolution  of  the  Data  Precision 
DVM  is  one  microvolt.  DVM  1  is  dedicated  to  reading  tne  thermocouple 
voltage  from  which  the  device  temperature  can  be  determined.  Einary- 
coded-decimal  data  from  DVM  1  and  DVM  2  is  fed  to  the  optically-isolated 
digital  input  cards  in  slots  407  through  41C. 

A  library  of  relocatable  subroutines  which  either  control  or 
receive  data  from  the  instruments  in  the  system  has  been  previously 
developed  by  the  author.  These  subroutines  make  the  job  of  writing 
data  acquisition  programs  much  less  tedious  and  reduce  the  chances  for 
program  errors. 

1.4  Temperature  Control  System  Hardware 

The  hardware  units  specifically  related  to  the  temperature  control 
system  are  the  D/A  voltage  converter  card,  the  HP  6224B  power  supply 
(Figure  1.1)  and  the  probe  which  holds  the  semiconductor  device-under¬ 
test  and  also  forms  the  thermal  portion  of  the  system  (Figure  1.2). 
Hereafter,  this  probe  will  be  referred  to  as  the  tnermal  plant.  The  D/A 
voltage  converter  is  capable  of  delivering  a  DC  voltage  in  the  range 
-10.240V  to  +10. 235V  in  5mV  steps.  It  is  used  to  voltage  program  the 
6224B  power  supply  with  a  gain  of  two.  Power  is  supplied  to  the  heater 
wire  of  the  thermal  plant  by  the  6224B.  Its  output  voltage  is  controlled 
by  the  software  subroutine  SETVH  which  also  limits  the  output  to  20V. 

A  listing  of  the  subroutine  SETVH  is  included  in  Appendix  A.  The  resis¬ 
tance  of  the  heater  wire  is  13*35  ohms.  Therefore,  the  maximum  power 
which  can  be  supplied  to  the  thermal  plant  is  aoout  30  Watts.  Since 
energy  cannot  be  removed  from  the  plant  via  the  heater  wire,  the  minimum 
power  is  zero  Watt3.  The  unipolarity  and  saturation  of  the  power  which 
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can  be  supplied  to  the  thermal  plant  result  in  a  limited  region  of  linear 
operation.  Small  quantization  errors  also  exist  in  the  power  supplied 
to  the  heater  and  the  measured  temperature  due  to  the  use  of  digital 
signals  to  represent  these  quantities. 

The  thermal  plant  was  designed  by  previous  members  of  the  Solid 
State  Electronics  Laboratory.  Cn  the  upper  copper  device  block  there 
is  a  socket  which  accepts  an  eight-pin  TO-5  can  which  contains  one  to 
four  solid  state  devices.  An  anodized  aluminum  cover  is  placed  over  the 
device  can  and  is  in  direct  thermal  contact  with  the  copper  block.  The 
layer  of  aluminum  oxide  prevents  any  electrical  contact  to  the  TO-5  can. 
The  measurement  thermocouple  junction  is  imbedded  in  the  aluminum  device 
cover  and  is  also  electrically  insulated.  Since  the  device  is  complete¬ 
ly  surrounded  by  the  thermal  mass  of  the  copper  block  and  the  aluminum 
cover,  the  temperature  measured  in  the  aluminum  cover  will  accurately 
reflect  the  temperature  of  the  semiconductor  device.  Due  to  this  large 
thermal  mass,  the  device  is  not  subject  to  rapid  temperature  fluctua¬ 
tions.  The  copper  block  at  the  lower  end  of  the  brass  rod  is  used  to 
reduce  the  sensitivity  of  the  plant  parameters  to  the  level  of  liquid 
nitrogen  in  the  dewar.  In  normal  use,  the  nitrogen  level  should  be 
even  with  the  top  of  the  lower  copper  block.  In  order  to  maintain  a 
constant  temperature  at  the  device  block,  the  power  supplied  through 
the  heater  wire  and  also  from  the  outside  environment  must  be  equal  to 
the  power  which  flows  out  of  the  device  block  through  the  brass  rod  and 
also  to  the  air  in  the  d°war. 

1.5  Temperature  Measurement  with  a  Coppe^-Constantan  Thermocouple 

When  two  dissimilar  metals  are  joined,  a  voltage  is  developed  across 
the  junction.  Thi3  phenomenon  is  called  the  Seebeck  effect  and  is  the 


8 


basis  of  the  operation  of  thermocouples.  A  thermocouple  consists  of  two 
wires  of  different  types  of  metals,  twisted  together,  and  heated  at  the 
junction  until  the  two  wires  begin  to  alloy.  A  Seebeck  voltage  is  gener¬ 
ated  at  the  junction  of  the  two  wires,  whose  magnitude  depends  on  the 
composition  of  the  metals  and  the  temperature  of  the  junction.  As  the 
temperature  of  the  junction  increases  so  does  the  voltage  at  the  open 
ends  of  the  wires.  Two  thermocouple  junctions  connected  in  series  are 
used  in  the  temperature  measuring  circuit.  When  one  of  the  junctions 
is  held  at  a  known,  constant  temperature,  the  voltage  produced  by  the 
pair  of  junctions  can  be  used  to  determine  the  temperature  of  the  non¬ 
reference  junction.  [31 

The  thermocouples  used  in  this  temperature  control  system  are 
type  T  (copper  vs.  constantan).  Constantan  is  a  copper-nickel  alloy. 

The  copper  wire  has  a  positive  potential  with  respect  to  the  constantan 
wire.  The  temperature  measurement  circuit  is  shown  in  Figure  1.3. 

Note  in  the  figure  that  junctions  of  dissimilar  metals  also  occur  where 
the  coaxial  cable  leading  to  the  input  of  the  digital  voltmeter  connects 
to  the  copper  thermocouple  wires.  These  junctions  are,  however,  in 
close  proximity  and  should  be  at  the  same  temperature  so  that  the  volt¬ 
ages  generated  at  these  junctions  will  cancel  each  other.  The  lower 
reference  junction  is  immersed  at  all  times  in  liquid  nitrogen  which  is 
at  its  boiling  point,  77. 3^8  K,  according  to  the  International  Practical 
Temperature  Scale  adopted  in  1968  (IPTS-68).  [4] 

The  type  T  thermocouple  is  the  only  one  of  the  standardized  types 
for  which  limits  of  error  are  set  below  0°C.  The  thermocouple  wires 
used  in  this  temperature  control  system  are  specified  by  the  manufac¬ 
turer  (Claud  S.  Gordon  Co.,  Richmond,  IL)  to  meet  the  special  limits  of 
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Figure  1.3  Temperature  measurement  circuit. 
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error  which  are  plus  or  minus  one  percent  of  the  measured  temperature 
between  -184  and  -59°C  and  plus  or  minus  0.8°C  between  -59  and  93°C 
when  the  reference  junction  is  at  0°C.  [4]  With  these  limits  in  mind, 
the  absolute  accuracy  of  the  measured  temperature  should  not  be  con¬ 
sidered  to  be  better  than  plus  or  minus  one  percent  of  the  measured 
temperature.  Calibration  of  the  thermocouple  against  a  platinum 
resistance  thermometer  may  improve  the  absolute  accuracy. 

A  fourteen  term  power  series  expansion  for  the  thermoelectric 
voltage  in  microvolts  of  type  T  thermocouples  is  given  in  Reference  [4], 
The  series  is  not  repeated  here.  This  power  series  with  temperature  in 
degrees  Celsius  as  the  independent  variable  is  for  the  reference  junc¬ 
tion  at  0°C.  To  adapt  this  series  to  the  reference  junction  at  the 
boiling  point  of  nitrogen,  all  that  need3  to  be  done  is  to  calculate  the 
thermoelectric  voltage  at  the  reference  temperature  (-195.802°C)  and 
subtract  this  value  from  the  voltage  obtained  from  the  series.  Of 
course,  the  temperature  used  as  the  argument  of  the  series  must  still  be 
in  degrees  Celsius.  The  temperature  may  be  determined  from  the  measured 
thermocouple  voltage  by  solving  the  equation 

f(T)  -  f(-195.802)  -  (thermocouple  voltage)  =  0  (1.1) 

where  f  denotes  the  power  series,  T  is  the  temperature  in  degrees 
Celsius,  and  the  thermocouple  voltage  is  in  microvolts.  The  function 
TEMP,  listed  with  the  control  program  in  Appendix  A,  uses  the  Newton- 
Raphson  Technique  to  solve  Equation  (1.1)  for  the  temperature  in  degrees 
Kelvin.  Function  TEMP  converges  in  less  than  seven  iterations. 

1.6  Approximate  Model  of  the  Thermal  Plant 

Before  one  can  go  about  designing  a  feedback  control  system,  one 
must  have  a  good  model  of  the  process  to  be  controlled.  In  order  to 
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formulate  a  model  of  the  thermal  plant,  step  response  data  was  taken 
for  several  input  power  levels.  The  mathematical  form  of  the  step 
responses  was  determined  using  a  fitting  program  based  on  a  subroutine 
called  CURFIT  from  Reference  [5].  The  data  was  successfully  fitted  to 
a  function  consisting  of  two  exponential  terms  and  a  constant  term  of 
the  form 


T(t)  =  A1exp(-A2t)  +  A^exp(-Ajjt)  +  .  (1.2) 

This  is  the  same  form  of  response  which  would  be  obtained  when  a  step 
function  is  the  input  to  a  system  with  transfer  function  cf  the  form 
T(s)  A(s  +  z1 ) 

G(s)  =  -  =  -  ,  (1.3) 

P(s)  (s  +  p 1 ) ( S  +  P2) 

where  T(s)  is  the  Laplace  transform  of  the  temperature  response  and  P(s) 
is  the  Laplace  transform  of  the  input  power.  From  the  information  pro¬ 
vided  by  the  fitted  parameters  of  Equation  (1.2),  the  following  approxi¬ 
mate  values  were  found  for  the  parameters  of  the  transfer  function  in 
Equation  (1.3):  A  =  0.015,  z ^  =  0.0002,  p.]  =  0.0012,  and  p^  r  0.00017. 

A  lumped  model  of  the  thermal  plant  was  obtained  from  this  approxi¬ 
mate  transfer  function.  Figure  1.4  shows  this  lumped  model  in 
electrical  analog  form.  The  values  of  all  the  elements  of  the  model 

except  R  were  obtained  by  equating  the  transfer  function  of  Equation 
o 

(1.3)  to  the  transfer  function  of  the  model,  given  by  Equation  (1.4) 
and  solving  for  the  parameters  of  the  model. 


T(s) 


Ph(s) 


(s  +  1/R  C  )/C. 

3  3  u 


(1.4) 


♦  s(1/R  C.  +  1/R  C 
r  d  a  a 


1/R  C.)  ♦  1/R  C  R  C. 
ad  a  a  r  d 


The  value  of  R  was  calculated  from  the  temperature  of  the  device  block 
o 


when  P. 


0  W  according  to  the  following. 


R  s  R  (T  -  T ) / ( T  -  77.348) 
o  r  o 


(1.5) 


t  ■ 

i 


130  K/VI 


Figure  1.4  Lumped  model  of  the  thermal  plant. 


When  a  brass  rod  is  used,  Rr  is  14.7  K/W  and  T  is  about  100  K  with 

P,  =  0.  If  T  is  taken  as  300  K  then  R  is  found  to  be  about  130  K/W. 
ho  o 

So  the  power  entering  the  system  from  the  outside  environment  is  given 
by 

P0  =  (Tq  -  T)/1 30  W.  (1.6) 

To  check  the  validity  of  this  model,  the  thermal  resistance  of  the 
brass  rod,  R^,  and  the  thermal  capacitance  of  the  device  block  and 
cover,  C^,  were  calculated  from  the  dimensions  and  the  thermal  conduc¬ 
tivity  or  specific  heat  capacity  of  these  elements.  The  thermal  resis¬ 
tance  of  the  brass  rod  assuming  no  heat  loss  from  the  sides  of  the  rod 
was  calculated  to  be  19.6  K/W.  Of  course  if  heat  loss  from  the  sides  of 
the  rod  was  considered,  the  calculated  value  would  be  significantly 
less.  Thus  the  experimental  value  of  14.7  K/W  for  Rr  in  the  model  is 
roughly  correct.  The  calculated  thermal  capacitance  of  the  device 
block  and  cover  assuming  no  temperature  gradients  within  the  block  was 
8C  J/K.  When  one  considers  that  temperature  gradients  do  exist  within 
the  device  block  during  the  rise  time  of  the  step  response,  it  is 
logical  that  the  apparent  thermal  capacitance  will  be  smaller  since  the 
outside  regions  of  the  device  block  where  the  temperature  is  measured 
are  heated  before  the  inside  of  the  block.  Less  power  will  be  required 
to  cause  a  given  change  in  temperature  in  a  given  time  near  the  outside 
of  the  block  than  if  the  temperature  were  measured  in  the  middle  of  the 
block.  In  non-equilibrium  conditions  the  effective  thermal  capacitance 
of  the  block  will  be  reduced.  When  this  point  is  considered,  the 
experimental  value  of  66.7  J/K  for  used  in  the  model  is  also  seen  to 
be  reasonable. 

It  must  be  pointed  out  that  this  model  is  only  approximate.  The 
use  of  a  lumped  model  for  a  system  which  is  obviously  distributed  is, 
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strictly  speaking,  incorrect.  There  is  a  time  delay  in  the  actual 
system  which  is  not  accounted  for  in  the  lumped  model.  However,  the 
simplification  of  the  design  procedures  which  results  from  using  a 
lumped  model  is  a  benefit  which  outweighs  the  disadvantage  of  the 
inaccuracy  in  the  model. 
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Chapter  2.  CONTROLLER  DESIGN  USING  DISCRETE-TIME  TECHNIQUES 

2. 1  Discrete-Time  Design  Procedure 

The  sampled  signals  present  within  the  temperature  control  system 
are  not  continuous  functions  of  time  and  are  difficult  to  handle  using 
Laplace  transforms.  The  z-transform  method  is  useful  in  simplifying  the 
analysis  and  design  of  sampled-data  systems.  This  technique  is  employed 
in  the  design  of  a  digital  discrete-time  controller  in  this  chapter. 

It  was  pointed  out  in  Chapter  1  that  the  temperature  control  system  has 
non-ideal  properties  such  as  quantization  errors  and  saturation  of  the 
power  which  can  be  supplied  to  the  thermal  plant.  These  characteristics 
will  be  neglected  in  order  to  simplify  the  analysis  and  design.  The 
quantization  errors  may  be  neglected  since  they  are  small  in  comparison 
to  the  magnitudes  of  the  signals  present  in  the  system.  The  saturation 
effect  in  the  power  available  to  the  heater  will  come  into  play  only  for 
large  temperature  steps  in  the  reference  input. 

The  design  procedures  described  below  follow  closely  those  given 
in  Reference  [6]  and  Chapter  9  of  Reference  [7].  Only  the  major  points 
used  in  arriving  at  a  design  algorithm  are  repeated  here. 

Consider  the  closed  loop  sampled-data  system  in  Figur?  2.1. 

The  transfer  function  of  the  closed  loop  system  is  given  b: 

C(z)  D(z)G(z) 

M(z)  =  -  =  -  (2.1) 

R(z)  1  +  D(z)G(z) 

If  M(z)  can  be  selected  so  that  the  closed  loop  system  meets  certain 
design  specifications  such  as  percent  overshoot  or  rise  time,  then-  the 
appropriate  digital  controller  transfer  function,  D(z),  can  be  found  by 

M(z)  1 


D(z) 


1  -  M(z) 


3(z) 


(2.2) 


synchronous  samplers 


R(z)  -  reference  input 
C(z)  -  controlled  variable 
E(z)  -  error  signal 

D(z)  -  digital  controller  transfer  function 
G(z)  -  plant  transfer  function 


Figure  2.1  Sampled-data  control  system  with  digital  controller 
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The  choice  of  M(z)  cannot  be  made  without  regard  to  the  transfer 
function  of  the  plant,  G(z).  The  controller  transfer  function  found 
from  Equation  (2.2)  must  be  physically  realizable  and  should  also  be 
stable.  Physical  realizability  of  D(z)  implies  that  the  degree  of  the 
numerator  must  not  exceed  that  of  the  denominator.  The  digital 
controller  will  be  stable  if  the  poles  of  D(z)  lie  inside  the  unit 
circle  in  the  z-plane. 

The  closed  loop  transfer  function,  M(z),  must  now  be  related  to  the 
performance  specifications  of  the  system.  Consider  M(z)  in  the  form 

K(z  -  z  ) 

M(z)  =  - !— — —  (2.3) 

(z  -  Pl)(z  -  p,) 

where  p1  and  fry  are  a  pair  of  complex  conjugate  poles.  K  is  chosen 
so  that  the  steady  state  error  for  a  step  input  is  zero.  By  application 
of  the  Final-Value  Theorem,  K  is  found  to  be 
(1  -  p  )  (1  -  p  ) 

K  =  - ! - —  (2.4) 

1  '  Z1 

The  step  response  of  the  closed  lcop  system  is  found  to  be 
c(nT)  =  1  +  seca  tp^l0  cos  (n$  +  a  -  r  )  for  n=0,1,2...  (2.5) 
where  b  =  arg(p1)  and  a  =  arg(p1  -  z  )  -  arg(p^  -  1)  +  V2 .  The 
complex  pole  p1  has  a  counterpart  in  the  s-plane  denoted  p'  such  that 
p1  =  exp(p'1T)  =  exp(-;^nT)  exp( Jai  T)  (2.6) 

where  i  is  the  damping  ratio  and  n>n  Is  the  natural  frequency  of  the 
system.  From  this  relation  it  is  seen  that 

Ip.,  I  =  exp(-s^'nT)  and  b  =  ^n*1  -  ;*  T  (2.7) 

When  these  relations  are  substituted  into  Equation  (2.5),  a  con¬ 
tinuous  function  which  passes  through  the  sampled  points  is  obtained  by 


reolacing  nT  with  t. 
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c ( t )  =  1  +  sec  a  exp(-ru  t)  cos(uj  /I  - 7 2  t  +  a  -•*)  (2.6) 

n  n 

This  equation  can  be  solved  for  the  time  at  which  the  maximum  response 
occurs . 


max 


UJ  /I 

n  s 


jarctan  f- 


-c 


/r 


(2.9) 


Substituting  this  expression  for  t  in  Equation  (2.8),  it  is  found 

max 


that  the  overshoot,  abbreviated  O.S.  ,  for  a  unit  step  input  is 


O.S.  =  sec  a  exp< - (arctan( — - — ■  |  -  a+rj }  (2. 10 ) 

/r '  /l 

For  a  given  damping  ratio,  Q,  the  value  of  a  which  gives  the 

minimum  overshoot  is  given  by 

a  =  arc  tan  (  -  <;/  /I  -  c,'^)  (2.11) 

If  a  is  chosen  by  (2.11),  using  (2.8)  Equation  (2.9)  may  be  written  as 

t  =  it  T  /  $  (2.12) 

max 

The  angle  of  the  complex  pole,  p^  may  be  found  from  Equation  (2.12) 

for  a  specified  t  .  The  magnitude  of  p  can  be  found  by  combining 
max  1 

the  two  relations  of  (2.7)  to  obtain 

Ip^  =  exp  (-i K  /  »'1  -  C2  )  (2.13) 

The  zero,  z^  lies  on  the  real  axis  of  the  z-plane  and  can  be 
found  from  the  relation 

a  =  arg(  p  -  z^  -  arg(p1  -  1)  ♦  -t/2  (2.14) 

The  constant  K  can  now  be  found  from  Equation  (2.4).  At  this 

point,  all  of  the  parameters  of  the  closed  loop  transfer  function 

have  been  found  for  any  choice  of  ■-  and  t  .  The  discrete-time 

max 

controller  transfer  function,  D(z),  can  be  found  from  Equation  (2.2) 
provided  the  transfer  function  of  the  plant,  G(z),  is  known.  The 
design  procedure  may  be  developed  into  a  point-by-point  design 
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algorithm  suitable  for  execution  by  a  computer  program.  This  is  j 

done  in  Appendix  B.  j 

2.2  Discrete-Time  Temperature  Controller  Design  ! 

The  design  procedures  developed  in  Section  2.1  and  Appendix  E 

) 

will  be  applied  to  the  temperature  control  system  in  this  section.  I 

I 

The  voltage  applied  to  the  heater  wire  of  the  thermal  plant  is  held 

constant  until  the  next  sampling  instant  when  a  new  voltage  is  set. 

Thus,  the  power  input  to  the  thermal  plant  is  subjected  to  a  sample- 

and-hold  operation  commonly  referred  to  as  the  zero-order  hold.  A 

block  diagram  portraying  this  action  is  shown  in  Figure  2.2  (a).  The 

z-transformed  counterpart  is  showr.  in  Figure  2.2  (b)  for  a  sampling 

period  of  two  seconds.  The  zero-order  hold  may  be  considered  as  a 

part  of  the  plant.  The  transfer  function  of  the  plant  with  the  zero- 

order  hold  included  is  denoted  by  G  (z). 

P 

With  this  knowledge  of  the  transfer  function  of  the  plant,  the 
design  algorithm  described  in  Appendix  B  may  be  applied  to  the  temper¬ 
ature  control  system  to  find  the  controller  transfer  function.  The 
controller  transfer  function  may  be  written  in  the  form 

P(z)  a  +  a  Z~  *  a  z"c  +  a,z"  j 

D  ( z )  =  ■  -  ~  "  g— ■  -- - 1  —  (2.15)  j 

E(z)  1  +  b.jZ  +  b^z  ♦  b-z 

where  P(z)  is  the  power  output  of  the  controller  and  E(z)  is  the  error 

-n 

signal  input  to  the  controller.  In  Eouation  (2.15),  the  factor  z 
implies  a  time  delay  of  nT  seconds.  D(z)  may  easily  be  implemented 
by  a  computer  program.  The  operation  of  the  computer  program  is 
depicted  in  block  diagram  form  in  Figure  2.?.  The  output  sequence  of 
the  digital  discrete-time  controller  depends  on  the  past  outputs  as 


well  as  the  oast  and  present  inputs. 


G  ( s ) 
P 


J 


(a) 


Figure  2.2  Block  diagram  of  thermal  plant  with  zero-order  saraple- 

and-hold.  (a)  Laplace  transform  diagram,  and  (b)  z-trans- 
form  diagram.  The  plant  transfer  function  was  found 
experimentally . 


Figure  2.3  Block  diagram  of  computer  program  implementation  of  D(z).  e  (t)  denotes  the  error 
Vernal  incut  and  p*(t)  denotes  the  power  output.  After  Kuo  [7]. 


Several  discrete-time  controllers  were  designed  with  tne  aid  of 
program  DSN,  listed  in  Appendix  B,  and  tested.  None  of  the  discrete¬ 
time  controllers  exhibited  totally  satisfactory  responses.  An  example 
of  the  step  response  of  one  of  the  discrete-time  controllers  is  shown  in 
Figure  2.4.  The  transfer  function  of  the  controller  which  produced  this 
response  is  given  in  the  sample  output  of  program  DSN  in  Appendix  B.  The 
response  shown  in  Figure  2.4  is  unsatisfactory  for  the  following  reasons: 

(i)  there  is  a  steady  state  error  of  0.5  K, 

(ii)  the  response  is  oscillatory  before  reaching  the  steady  state, 

(iii)  the  settling  time  of  about  30  minutes  is  excessive. 

The  most  serious  of  these  deficiencies  is  the  non-zero  steady  state 
error.  A  deeper  consideration  of  the  design  algorithm  reveals  why  the 
steady  state  error  is  not  zero.  Referring  to  Equation  (2.2)  it  is  seen 
that  the  determination  of  D(z)  which  shoulc  give  the  desired  closed  lccp 
transfer  function,  M(z),  depends  on  a  rather  precise  knowledge  of  the 
plant  transfer  function,  G(z).  In  the  case  of  the  temperature  control 
system,  G(z)  is  only  approximately  known.  Ihe  constant  K  of  the  closed 
loop  transfer  function  was  chosen  for  zero  steady  state  error  using 
Equation  (2.4).  The  values  used  for  the  poles  and  zero  in  this 
equation  are  calculated  values  which  may  not  be  equal  to  the  actual 
poles  or  zero  of  the  closed  loop  system  due  to  the  inexact  knowledge 
of  the  plant  transfer  function.  Therefore,  D(z)  found  from  Equation 
(2.2)  may  not  yield  the  desired  closed  loop  transfer  function  which 
would  produce  a  response  with  zero  steady  state  error.  Since  this 
particular  discrete-time  design  approach  coes  not  produce  a  closed 
Loop  system  which  meets  the  design  goals,  some  other  techniques  must 
be  used  to  design  the  control  system.  In  Chapter  3  the  control  system 
is  treated  as  an  approximation  to  a  continuous-time  system. 
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TEMPERATURE  CONTROLLER  STEP  RESPONSE 


CONTROLLER  VERSION  CN06 
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Figure  2.4  Response  of  temperature  control  system  with  a  digital  controller 
to  a  step  input  from  150  to  155  K. 


«  • 
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Chapter  3.  CONTROLLER  DESIGN  USING  CONTINUOUS-TIME  TECHNIQUES 

3.1  Proportional.  Integral,  and  Proportional  plus  Integral  Control 

Consider  a  closed  loop  system  of  the  form  shown  in  Figure  3.1(a). 

If  the  sampling  period,  T  ,  is  small  compared  to  the  system  time  con- 

3 

stants,  the  system  may  be  treated  as  a  continuous-time  system  for  the 
purposes  of  analysis  and  design  (Figure  3.1(b)).  This  approach  is 
applied  to  the  temperature  control  system  in  this  chapter. 

The  plant  transfer  function  was  found  in  Chapter  l.  The  controller 
transfer  function,  D(s),  must  be  determined  so  that  the  desired  response 
of  the  closed  loop  system  can  be  obtained.  In  this  temperature  control 
system,  D(s)  will  be  implemented  using  a  digital  computer  (HP  1000M) . 

The  saturation  effect  in  the  power  which  can  be  supplied  to  the  thermal 
plant  is  neglected  in  the  initial  design.  The  design  goals  are: 

(1)  the  system  should  be  stable, 

(2)  the  steady  state  error  for  a  step  input  must  be  zero,  and 

(3)  the  step  response  overshoot  and  settling  time  should  be 
reasonable . 

The  question  which  arises  now  is  what  form  the  controller  transfer 
function  should  take  so  that  the  design  goals  are  met. 

The  plant  transfer  function  is  of  the  form 

A(s  ♦  z.) 

G(s)  =  - - -  (3.1) 

(s  +  p^Hs  +  p2) 

where  A  =  0.015,  z 1  =  C.0002,  p1  =  0.0012,  and  p2  =  0.00017. 

If  proportional  control  is  used  the  controller  transfer  function  will 
simply  be  a  fixed  gain,  Kp,  as  shown  in  Figure  3.2(3).  The  closed  loop 
transfer  function  is  of  the  form 


Figure  3.1  (a)  Discrete-time  representation  of  the  temperature 

control  system,  and  (b)  Continuous-time  representation 
of  the  temperature  control  system. 


(a) 


(c) 

Figure  3.2  Temperature  control  system  block  diagram  with 
(a)  proportional  control,  (b)  integral  control, 
(c)  proportional  plus  integral  control. 
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T(s) 


AKp(s  z.,) 


(3-2) 


P(s)  s  +  (p1  +  p2  ♦  AKp)s  +  p^p  *  AKp21 
The  Routh-Hurwitz  stability  test  shows  that  this  system  will  be  stable 
for  Kp  >  0.  The  steady  state  error  for  a  unit  step  input  i3  non-zero 
for  finite  values  of  Kp.  [3] 

Egg  =  P-jP2  >  (P,P2  ♦  AKpZ-j )  (3.3) 

It  is  seen  that  the  proportional  controller  cannot  meet  design  goal 
number  2. 

If  integral  control  is  used  the  controller  transfer  function  will 
take  the  form  K^/s  as  shown  in  Figure  3.2(b).  The  closec  loop  transfer 
function  becomes 


T(s) 


AK  (a  +  2 ^ ) 


3 2 

R(s)  s  «■  (p1  +  p_,)s  ♦  (P1P2  +  AKj)s  ♦  AKjZ 


(3.4) 


This  system  will  be  stable  for  >  0  according  to  the  Routh-Hurwitz 
test.  The  steady  state  error  for  a  unit  step  input  is  found  to  be  zero 
by  application  of  the  Final-Value  Theorem.  Integral  control  was  tested 
on  the  control  system  and  it  was  found  that  design  voai  number  3  could 
not  be  satisfied.  The  overshoot  of  the  step  response  was  excessive  and 
the  settling  time  was  much  too  long. 

If  proportional  plus  integral  control  are  employed,  the  controller 
transfer  function  will  be  of  the  form  shown  in  Figure  3.2(c).  The 
closed  loop  transfer  function  is 

T(s) 


AKp(s  ♦  KT/Kp)(s  +  z ^ ) 


P(s)  +  (pi  ♦  p^  AKp)s2  ( p 1 +  AK^  *  AKpZ  ^ ) s  ♦  AK^ 


(3.5) 


For  p^  p2  -  z1  >  0,  the  system  is  stable  for  Kp  >  0  and  >  0.  These 
conditions  are  met  by  the  temperature  control  system.  The  steady  state 
error  for  a  unit  step  input  is  zero.  Thus,  a  proportional  plus  integral 


J 


controller  meets  design  goals  1  and  2.  Root  locus  and  Laplace  transform 


techniques  can  be  used  to  study  the  form  of  responses  which  will  be 
obtained  as  Kp  and  are  varied.  Proportional  plus  integral  control  is 
definitely  the  best  choice  for  the  temperature  control  system. 

3.2  Design  Using  Root  Locus  Techniques 

Information  on  how  the  proportional  constant,  Kp,  and  the  integral 
constant,  K^,  should  be  chosen  may  be  obtained  from  the  loci  of  the 
poles  of  the  closed  loop  transfer  function.  To  get  this  information 
Kp  can  be  fixed  at  a  chosen  value  and  varied  from  zero  upward.  The 
characteristic  equation  is  obtained  by  setting  the  denominator  of 
Equation  (3-5)  equal  to  zero. 

s^  +  (p1  +  p2  +  AKp)s2  +  (p1P2  +  AKj  +  AKpZ^s  +  AKjz-)  ~  0  (3.6) 

The  zeros  of  the  characteristic  equation  wnich  are  the  poles  of  the 
transfer  function  are  found  and  plotted  in  the  s-plane  for  each  value 
of  Kp  and  K^.. 

Figure  3-3  shows  the  behavior  of  the  poles  in  the  s-plane  for 
Kp  =  2.0  and  K  varied  from  0  to  0.05  in  steps  of  0.0001.  The  behavior 
of  the  poles  near  the  closed  loop  zero  at  s  =  -0.0002  is  shown  in  Figure 
3.4.  This  feature  is  not  visible  in  Figure  3-3  due  to  the  plotting 
scale.  The  points  on  the  real  axis  where  two  poles  meet  and  become  a 
pair  of  complex  conjugate  poles  are  called  breakaway  points.  As  the 
fixed  value  of  Kp  is  increased,  the  breakaway  point  shown  in  Figure  3.5 
moves  to  the  left  in  the  s-plane.  One  of  the  poles  approaches  the  zero 
at  s  s  -0.0002  very  rapidly  as  is  varied  so  that  there  will  be  a 
pole-zero  cancellation  except  when  is  very  small.  For  Kp  =  2.0,  the 

breakaway  point  occurs  at  s  =  -0.0156  when  s  0.0162.  For  K ^  <  0.0162, 

the  two  uncancelled  poles  will  be  on  the  real  axis  and  the  step  response 
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Figure  3.3  Loci  of  poles  of  temperature  control  system  with  a 
proportional  plus  integral  controller.  K  =  2.0, 
.<T  varied  from  0.00  to  0.05  in  steps  of  0.0001. 
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Behavior  of  poles  near  the  zero  at  s  :  -0.0002. 
Kp  =  2.0,  Kt  varied  in  steps  of  0.000002. 
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will  consist  cf  two  decaying  exponential  terms  plus  a  constant.  For 
example,  when  Kp  =  2.0  and  Kj  =  0.01  the  response  to  a  unit  steD  input 
is 

T(t)  =  1  -  1.248  exp(-0. 02522  t)  +  0.24?  exp( -0 . 005945  t)  (3.7) 

For  Xj  >  0.0162  the  step  response  will  consist  of  a  damped  sinusoid 
plus  a  constant  term. 

The  methods  described  above  can  be  used  to  investigate  the 
theoretical  system  responses  as  Kp  and  are  varied,  but  in  order  to 
choose  the  values  of  Kp  and  which  yield  the  best  performance  for  the 
actual  system,  several  controllers  must  be  implemented  and  tested  with 
the  system. 

3.3  Computer  Implementation  of  the  Temperature  Controller 

The  transfer  function  of  the  proportional  plus  integral  controller 
may  be  implemented  using  a  computer  program.  The  reference  input  to  the 
controller  is  passed  to  it  in  the  form  of  digital  data  by  a  monitor  cr 
or  data  acquisition  program  by  means  of  Class  1/0  calls  which  permit 
communication  between  separate  programs.  Class  I/O  is  provided  in  the 
RTE-IV  operating  system.  [9]  Class  1/0  communication  between  the  two 
programs  i3  illustrated  in  Figure  3-5. 

The  sampling  period  is  also  passed  to  the  controller  program  along 
with  the  reference  temperature.  This  flexibility  in  choosing  the 
sampling  period  is  a  definite  advantage.  A  data  acquisition  program 
might  instruct  the  temperature  control  program  to  run  every  two  seconds. 
When  the  data  acquisition  program  reached  a  point  where  the  sampling  of 
data  on  the  semiconductor  device  became  very  time-critical,  it  could 
set  the  sampling  period  of  the  temperature  control  program  to  some 


larger  value  so  that  the  temperature  control  program  would  not  use  any 
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CPU  time  during  this  critical  period.  After  the  data  has  been  taken, 
the  temperature  sampling  period  could  be  reset  to  two  seconds. 

In  the  temperature  control  program,  the  trapezoidal  approximaticn 
is  used  to  compute  the  integral  of  the  error  signal.  A  flow  chart 
showing  the  operations  performed  by  the  temperature  control  program 
is  shown  in  Figure  3.6.  The  temperature  control  program  is  listed  in 
Appendix  A  along  with  a  sample  monitor  program  which  serves  as  an 
example  of  how  to  use  the  control  program. 

3.9  Testing  of  the  Temperature  Control  System 

The  actual  response  of  the  temperature  control  system  to  a  step 

input  from  195  K  to  200  K,  with  K  =  2.0  and  KT  =  0.01  is  shown  in 

Pi 

Figure  3*7.  The  theoretical  step  response  for  the  same  values  of 
and  Kj  and  the  same  input  is  shown  in  Figure  3.8.  The  overshoot  of  the 
actual  response  is  due  to  the  time  delay  in  the  thermal  Diant  which  is 
not  accounted  for  in  the  lumped  model.  The  best  values  of  Kp  anc  K. 
were  found  to  be  K  =  2.0  and  <(T  =  0.01  with  a  brass  rod  used  in  the 

r  x 

thermal  plant. 

It  was  mentioned  in  Chapter  1  chat  the  power  which  can  be  supplied 
to  the  thermal  plant  saturates  at  about  30  VI.  For  steps  in  the  input 
larger  than  15  K,  saturation  of  the  applied  power  will  occur.  Large 
overshoots  and  long  settling  times  are  the  results  of  this  nonlinearity. 
A  simple  correction  can  be  made  to  compensate  for  this  nonlinearity. 
Whenever  the  magnitude  of  the  error  signal  is  greater  than  15  K,  the 
integral  of  the  error  signal  is  not  updated.  Response  to  a  53  K  step  m 
the  input  without  this  correction  is  shown  in  Figure  3.9.  The  improve¬ 
ment  of  the  response  with  the  corr°ction  for  the  nonlinearity  included 
is  very  significant  as  shown  in  Figure  j.  1C. 


perform  a  Class  1/0  Get  to  receive  data 
from  monitor  or  data  acquisition  program 


content  of  A-register  >  0 


e  s 


retrieve  data  from  buffer  (tref,  deltat) 


compute  a  new  delay  factor  from  deltat 


read  DVM  to  get  thermocouple  voltage 


find  temperature  from  thermocouple  voltage 


compute  error  signal 


A^s7~eTror  s  ^gTaT)~~^~l*5  T 


compute  integral  of  error  signal 


compute  a  new  power 


i— °  - ■ 


voltage  =  0 


voltage  =  y/power*resistance 


set  heater  voltage 


wait  for 

delay  time 

Fi  ure  3.6  Flow  chart  of  emperature  control  program 
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THEORETICAL  step  response 
CONTROLLER  VERSION  CN22 


0  600  1200  1800  2400  3000  3600 

TIME  (S) 


figure  3.o  Theoretical  step  response  based  on  tne  model  of  the  thermal 
plant  for  the  same  conditions  as  Figure  3.7. 
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Figure  3-9  Response  to  a  50  K  step  input  without  the  correction  for 
saturation  of  heater  power. 
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Fiaure  3. 10  Response  to  3  50  K  step  input  with  the  correction  for 
saturation  of  heater  power. 
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The  temperature  control  system  gives  very  satisfactory  performance 


over  the  temperature  range  from  100  K  to  300  K  with  a  brass  rod  used 
between  the  device  block  and  the  lower  copper  block  of  the  thermal  plant 
The  temperature  in  this  range  can  be  controlled  with  a  precision  of 
better  than  plus  or  minus  0.1  K.  In  order  to  obtain  stable  temperatures 
lower  than  100  K,  a  copper  rod  which  has  a  smaller  thermal  resistance 
must  be  substituted  for  the  brass  rod.  The  copper  rod  has  a  length  of 
7.635  cm.  and  a  diameter  of  1.275  cm.  The  thermal  resistance  of  tnis 
rod  is  calculated  to  be  1.46  K/W.  When  this  thermal  resistance  is 
substituted  into  the  model  and  the  design  is  performed  and  tested  as 

described  earlier,  it  is  found  that  Kp  =  2.0  and  K  =  0.01  gives  very 

good  results.  The  theory  shows  that  0.02  would  be  a  better  value  for  KT 
However,  experiments  showed  that  K^.  =  0.01  gave  a  better  steady  state 
response.  Accumulation  of  the  integral  of  the  error  signal  has  an 
averaging  effect  on  any  noise  in  the  system.  As  the  integral  control 
constant  is  increased,  this  averaging  effect  is  diminished  and  response 
to  the  noise  in  the  system  becomes  apparent  in  the  form  of  small  oscil¬ 
lations  in  the  temperature.  is  chosen  equal  to  0.01  to  obtain  better 

immunity  to  noi3e.  The  minimum  stable  temperature  which  can  be  main¬ 
tained  is  about  85  K.  The  precision  of  the  temperature  control  is  the 
same  as  that  for  the  brass  rod. 

Of  course  there  is  a  limit  on  how  long  a  certain  temperature  can 
be  maintained  with  either  type  of  rod  due  to  the  evaporation  of  the 
liquid  nitrogen.  Usually  the  liquid  nitrogen  will  last  six  to  eight 
hours  before  it  must  be  replenished. 

Several  sampling  periods  were  tested.  The  temperature  control 
system  performed  well  with  sampling  periods  less  than  or  equal  to 


»'  • 
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fifteen  seconds. 


AO 

The  response  to  a  ramp  input  R(t)  =  100  0.05t  K  where  t  is  in 

seconds  is  shown  in  Figure  3. 11.  The  ramp  response  is  of  interest  if 
Thermally  Stimulated  Capacitance  (TSCAP)  measurements  are  to  be  taken.  [1] 
The  response  to  faster  ramp  inputs  is  nonlinear  due  to  the  saturation 
of  the  power  supplied  to  the  heater  wire  of  the  thermal  plant. 


TEMPERRTURE  (K) 

100  130  160  190  220  250  280 
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TEMPERATURE  CONTROLLER  RAMP  RESPONSE 
CONTROLLER  VERSION  CN22 


Figure  3-11  Response  to  a  ramp  input,  R(t)  a  100  +  0.05t  K. 
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Chapter  4.  CONCLUSION 


The  computer-automated  temperature  control  system  described  in 
this  thesis  is  a  significant  step  toward  complete  automation  of  data 
acquisition  for  the  characterization  of  fundamental  electrical  proper¬ 
ties  of  semiconductor  materials. 

The  completion  of  this  project  required  work  in  a  wide  range  of 
topics  including  system  identification  and  modeling,  control  system 
design  using  both  discrete  and  continuous-time  techniques,  and  computer 
implementation  of  the  temperature  controller. 

A  design  algorithm  was  developed  for  a  controller  using  discrete¬ 
time  techniques.  This  algorithm  did  not  produce  a  controller  which 
performed  satisfactorily  because  it  required  knowledge  of  the  plant 
transfer  function  with  more  accuracy  than  was  obtainable.  A  controller 
which  met  all  of  the  system  design  goals  was  designed  by  application  of 
continuous-time  techniques.  This  controller  used  both  proportional  and 
integral  control.  Strictly  speaking,  the  temperature  control  system  is 
a  discrete-time  system.  However,  since  the  sampling  period  is  very 
small  comparec  to  the  system  time  constants,  this  system  can  be  approx¬ 
imated  by  a  cc ntinuous-tirae  system. 

The  final  implementation  of  the  temperature  control  system  is 
capable  of  controlling  temperature,  with  a  precision  of  +0.1  K,  in  the 
range  from  100  K  to  300  K  using  a  brass  rod  in  the  thermal  plant.  A 
copper  rod  must  be  substituted  for  the  brass  rod  to  obtain  temperatures 
below  100  K.  The  temperature  may  be  maintained  within  *0.1  .<  of  the 
reference  input  for  temperatures  down  to  85  K  using  the  copper  rod. 

The  absolute  accuracy  of  the  measured  temperature  is  determined  by  the 
manufacturer's  specifications  for  the  thermocouple  wire,  which  is  plus 


or  minus  one  percent  of  the  measured  temperature.  Better  absolute 
accuracy  of  the  measured  temperature  may  be  obtained  by  calibrating  the 
thermocouple  against  a  platinum  resistor  which  is  a  very  stable  and 
linear  temperature  transducer.  This  method  has  been  used  by  members  of 
our  laboratory  to  obtain  absolute  accuracy  of  ±0. 1  X. 
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APPENDIX  A. 


Listing  of  Temperature  Control  Program  and  Monitor  Program 
FTN4.L 

c 

C**“»  RNS790302  -  51  SOURCE  FILE  :  RNTCON 
C 

n  lltftlllllliftllft 

c  *  ■ 

PROGRAM  TCON 

C  »  * 

C 

THIS  IS  THE  TEMPERATURE  CONTROL  PROGRAM  TO  BE  USED 
C  WITH  A  BRASS  OR  COPPER  ROD  IN  THE  THERMAL  PLANT  CF  THE 

C  SYSTEM.  THIS  PROGRAM  IS  FOR  USE  ON  THE  MULTIPROGRAMMER  69*»0B 

C  MAINFRAME  (UNIT  ADDRESS  00). 

C 

C“M*  THIS  PROGRAM  COMMUNICATES  WITH  A  MONITOR  OR  DATA 
C  ACQUISITION  PROGRAM  VIA  CLASS  I/O  CALLS  PROVIDED  BY 

C  THE  RTE  IV  OPERATING  SYSTEM.  THIS  PROGRAM  RECEIVES 

C  THE  REFERENCE  TEMPERATURE,  TREF,  AND  THE  SAMPLING 

C  PERIOD,  DELTAT,  FROM  THE  MONITOR  OR  DATA  ACQUISITION 

C  PROGRAM.  THE  MAXIMUM  SAMPLING  PERIOD  SHOULD  BE  15  SEC. 

C  THE  MONITOR  OR  DATA  ACQUISITION  PROGRAM  HAS  THE  JOB  CF 

OBTAINING  A  CLASS  NUMBER  AND  SCHEDULING  THE  TEMPERATURE 
CONTROL  PROGRAM.  THE  MONITOR  OR  DATA  ACQUISITION  PROGRAM 
C  MUST  ALSO  RELEASE  THE  CLASS  NUMBER  AND  KILL  THE  CONTROL 

C  PROGRAM  WHEN  IT  IS  DESIRED  TO  END  THE  TEMPERATURE  CONTROL. 

C  THIS  ROUTINE  TAKES  APPROXIMATELY  JO  MILLISECONDS  TO  RUN. 

C 

C***»*  IMPORTANT  VARIABLES 


C  TREF  —  REFERENCE  TEMPERATURE  IN  DEGREES  KELVIN 

C  TMEAS  —  MEASURED  TEMPERATURE  IN  DEGREES  KELVIN 

C  DELTAT  -  SAMPLING  PERIOD 

C  ERRSIG  -  ERROR  SIGNAL  (TREF  -  TMEAS) 

C  CLDERR  -  LAST  VALUE  OF  ERRSIG 

C  INTERR  -  INTEGRAL  OF  THE  ERROR  SIGNAL 

C  TCV  -  THERMOCOUPLE  VOLTAGE 

C  KP - PROPORTIONAL  GAIN 

C  XI - INTEGRAL  GAIN 

C  RH  -  HEATER  RESISTANCE 

C  VH  -  VOLTAGE  APPLIED  TO  HEATER 

C  POWER  —  POWER  TO  BE  APPLIED  TO  THERMAL  PLANT 

C 

INTEGER  CLASS,  A,  B,  DELAY,  ERR,  PA  RAM ( 5 ) 


REAL  TREF,  DELTAT,  TMEAS,  ERRSIG,  GLDERR ,  INTERR 
REAL  TCV,  VH,  RH,  POWER,  KP ,  KI ,  MESSAG (2) 

DATA  KP/2.0/,  KI/1.0E-2/,  RH/ 1 3 -  35/ 

r 

I .  GET  THE  CLASS  NUMBER 

CALL  RMPAR ( PARAM ) 


f- 
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CLASS  =  PARAMO) 

C 

SET  BITS  15  AND  13  IN  CLASS  WORD 
C  (NO  WAIT  EIT  AND  DE-ALLOCATE  BIT) 

CLASS  =  IOR  (CLASS,  120000B) 

n 

INITIALIZE  OLDERR  AND  INTERR 
OLDERR  =0.0 
INTERR  =0.0 

C 

do  ft  CLASS  I/O  GET  TO  RETRIEVE  TREE  AND  DELTAT 
100  CALL  EXEC (21,  CLASS,  MESSAG,  4) 

C 

C***»«  CHECK  CONTENT  OF  A-REGISTER  TO  SEE  IF  NEW  DATA 
C  WAS  RECEIVED 

CALL  ABREG  (A,  B) 

IF  (A  -LT.  0)  GOTO  200 
C 

C»****  NEW  DATA  WAS  OBTAINED  IF  A  >=  0 
TREF  =  MESSAG ( 1 ) 

DELTAT  =  MESSAG ( 2  ) 

r 

COMPUTE  THE  DELAY  FROM  DELTAT 
DELAY  =  IFIX  (100. 0*DELTAT  -  4.0) 

C 

£#*«**  READ  the  TEMPERATURE  AND  COMPUTE  ERRSIG 
200  CALL  RDVM  (TCV) 

TMEAS  =  TEMP  (TCV) 

ERRSIG  =  TREF  -  TMEAS 
C 

C**»*»  CORRECTION  FOR  NONLINEARITY 

IF  (AES(ERRSIG)  .GT.  15.0)  GOTO  500 
C 

C#**#*  UPDATE  THE  INTEGRAL  OF  THE  ERROR  SIGNAL 

INTERR  =  INTERR  *■  (ERRSIG  +  0LDERR)/2.0  *  DELTAT 
C 

COMPUTE  THE  NEW  POWER  AND  SET  HEATER  VOLTAGE 
500  POWER  =  KP  *  ERRSIG  +  KI  *  INTERR 
C 

<;•••«•  powER  IS  UNIPOLAR 

IF  (POWER  .LT.  C.O)  GOTO  10 
VH  =  SORT  (POWER  *  RH) 

GOTO  20 

10  VH  =  0.0 

20  CALL  SETVH  (VH) 

r 

:»*»*»  SAVE  VALUE  OF  ERRSIG 
OLDERR  =  ERRSIG 
C 

C»*»**  WAIT  FOR  DELAY»1C  MILLISECONDS 
CALL  WAIT  (DELAY,  0,  ERR) 

GOTO  100 
END 


oooo  ooorsooonn  nnnoono 
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C _ * _ • _ # _ # _ * _ i _ • _ • _ • _ » _ • _ » _ • _ * _ t _ * _ 

C 

SUBROUTINE  RDVM  (TV) 

C 

C - PURPOSE - 

THIS  SUBROUTINE  READS  TWO  12  BIT  ECD  ENCODED  WORDS  FROM  THE 
DVM  THROUGH  THE  ISOLATED  INPUT  CARDS  IN  SLOTS  407  AND  408. 
THESE  TWO  DATA  WORDS  ARE  THEN  DECODED  TO  OBTAIN  THE  VOLTAGE 
READING. 

PARAMETERS - 

TV— STANDS  FOR  THERMOCOUPLE  VOLTAGE:  RETURNED  TO  MAIN  PRGM . 
DIMENSION  IBUFR(7 ) 

EUFFER  ELEMENTS 

17024GB— CONTROL  WORD:  ISL,SYE  ON 

1— PRIORITY  INDEX:  1  MEANS  START  AT  TCP  OF  LIST 
70000B— ADDRESS  SLOT  7 
100000B— ADDRESS  SLOT  3 
177777B— STOPPER 

IBUFR(6)  AND  IEUFR(?>  ARE  FOR  READ  IN  AREA. 

DATA  IEUFR/ 1 70240B , 1 , 70000B , 100000B, 1777778/ 

MAKE  A  POLL  ALL  INPUTS  REQUEST. 

ICCJE=1  CONTROL  WORDs 124B  EUFFER  LENGTH=7  IFUNC=5 

CALL  EXEC  ( 1 , 124B,IBUFR,7,5) 

C 

C  THE  REST  OF  THE  SUBROUTINE  DECODES  THE  TWO  BCD  ENCODED 

C  WORDS  TO  GET  THE  CORRECT  VOLTAGE  READING. 

C  FIRST,  OBTAIN  SIX  INTEGER  DIGITS,  THEN  THE  SIGN  INDICATOR 

C  (ISN)  AND  SCALE  INDICATOR  (IX). 

C 

IB0=IAND( i 7E , IBUFR(6 ) ) 

IB  1 =IAND(360B, IBUFR(6 ) )/l6 

IB2=IAND(7400B,IEUFR(6))/256 

IB3=IAND(17B,IBUFR(7)) 

IE4=IAND(360E,IBUFR(7))/16 
IB5=IAND(400B,IBUFR(7))/256 
ISN=IAND( W008,1BUFR(7)) 

IX=IAND(6000B,IBUFR(7))/1024-1 
IF(ISN.EQ.O)  GOTO  10 
SGN  =  1 . 

GOTO  20 
10  SGN=-  1 . 

C  NEXT,  COMBINE  INTEGER  DIGITS  AND  SIGN  TO  FORM  A  REAL  NUMBER. 

C  KEEPING  ALL  TERMS> 1  AND  USING  A  NESTED  POLYNOMIAL  HELPS 

C  AVOID  TRUNCATION  ERRORS. 

C 

20  :v=(  (  ( ( (IE5*1  .EUIB4)»1  .EUIE3)*1  .E1+IB2)»1  .EUIE1 ) *  1 .  E 1+IBO )  *SGN 


MAKE  SCALE  CORRECTION. 


1 
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IF  CX)  1,2,2 
C  .IV  SCALE 

1  TV=rV/1.0E5 
GOTO  c 

C  1  V  SCALE 

2  7V=rV/1.0E5 
GCT3  5 

3  IF  (IX. EC. 2)  GOTO  4 

C  1C  V  SCALE 

TV=TV/1 .0£4 
GOTO  5 

C  ICG  V  SCALE 

4  TVrTV/I.OE? 

5  RETURN 
END 

Q _ # _ » _ • _ * _ • _ * _ * _ * _ * _ * _ 

C 

FUNCTION  TEMP(TCV) 

C  VERSION  2.0  7/27/79  SEE  NOTEBOOK  1 6 1  P.  98 

C  *»•**  THIS  FUNCTION  FINDS  THE  TEMPERATURE  IN  KELVINS  FROM  THE  THERMC- 


C  »••**  COUPLE  VOLTAGE  REFERENCED  TO  THE  TEMPERATURE  OF  LIQUID  NITROGEN. 

C  •«***  THE  FUNCTION  USES  NEWTON'S  METHOD  TO  FIND  THE  ROOT  OF  THE  EQUATION 


C  #**•*  OF  THE  FORM  :  POLYNOMIAL  IN  T  -  THERMOCOUPLE  VOLTAGE  =  0, 

C  •*»**  WHERE  THE  POLYNOMIAL  IN  T  IS  A  14  TERM  POLYNOMIAL  BASED  ON  IPTS-68 
C  *****  DATA  FOR  TYPE  T  ( COPPER -CCNSTANTAN )  THERMOCOUPLES. 

DOUBLE  PRECISION  A , B , C , D , E , F , G , H , I , J , K , L , M , N , TC , TOLD , TNEW , EMFK 
DOUBLE  PRECISION  DSDT, ELN2 , EMF , C2 , C? , C4 , 05 , 06 , 07 , C8 , C9 , 01 0 
DOUBLE  PRECISION  C i  1 , C 12 , C 1 2 , C  1 4 
INTEGER  COUNT 

DATA  A / 3. 874e77?340D1/, 3/4. 4i 23932*820-2/, C/1. 140523S4g3D-4/ 

DATA  D/1 . 9974<*065cSD-5/, E/9. 0445401 1 S7D-7/ , F/2 . 27660 1 S5D4D-S/ 

DATA  G/2. 6247 4093800- 1 0/, H/3. 864292420 ID- 12/, 1/2. 62986 785 19D- 14/ 

DATA  J/i . 4261 3232 4od- l 6/, K/4. 883 3254 364D- 1 9/, L/ 1.08034746830-21/ 

DATA  M / 1 . 394Q29 1026D-24/ ,N/7 . 9795693 1 560-26/ , ELN2/0 . 5535590250D4/ 
DATA  C2/2 . ODO, ,C?/; . 3D0/ ,CU/4 . ODO/ ,C5/5 .ODO/ ,C6/6 . ODO/ , C7/7 . CDO/ 

DATA  C6/3.CDC7 .C9/9.0D0/ ,C 10/10. ODD/, Cl  1/1 1 . ODD/ , C 1 2/ 1 2 . ODO/ 

DATA  Cl  •/'.'.ODO/, 014/14. 000/ 

C 

c  CONVERT  FROM  UNITS  OF  VOLTS  TO  UNITS  OF  MICROVOLTS 

EMF =DBLE ( TCV ) *  1 .  )C6 
TNEW  =  1 . 0D2 

IF  vEMr .ST. 3.003)  TNEW=3-0D2 
DO  190  COUNT  s 1 . 1 ? 

TCLDsTNEW 
TCsTNEW-273. 1800 

SMFKs((  C  C  C  '  N*TC*m:  *T>L)*TC+K)*TC*J)*TC«-I  )*TC+.H)*TC+G)*TC*F 

4;  *TC+E  )*TC  -D  1  *TC+G  :  *TC*B,'  *TC  +  A  )  *TC*ELN2-EMF 
DEDT  =  (  ( (  (  '.  I  ( '  •  ;-4«N«TC*C1  ••m:  «TC*:12*L)*TC*C1-  »N)*TC*C10*J)*TC+ 

iC5»:  ) »TC*C6«H )  «TC*C7»3  )*TC-*-C6»F)  •TC*C“«E  )*TC^C4*D  ;  *TC+C2*C) *TC+C2* 

4F ;*TC+A 


o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o 
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TNSW=TCLD-EMFK/DEDT 
IF(DABS(TNEW-TOLD).LE.1.0D-5)  GOTO  1 
100  CONTINUE 

WRITE  ( 1 , 10) 

10  FORMAT  ("  FUNCTION  TEMP  DID  NOT  CONVERGE  TEMP  SET  TO  0.0  ") 
TNEW=0 . 0D0 

1  TEMP=SNGL(TNEW) 

RETURN 

END 

r* 

U 

Q _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ ft _ 

C 

SUBROUTINE  SETVH  (VH) 

C 

C  *»**»  RNS790329  -  U  SOURCE  FILE:  &MPLIB  BINARY  FILE:  JMPLIB 
PURPOSE - 

THIS  SUBROUTINE  SETS  THE  VOLTAGE  APPLIED  TO  THE  HEATER  COIL 
OF  THE  TEMPERATURE  CONTROLLER  FROM  THE  HP  6224B  POWER  SUPPLY 
BY  SETTING  THE  PROGRAMMING  VOLTAGE  GENERATED  BY  THE  A/D  CARD 
IN  SLOT  2  OF  THE  MULTIPROGRAMMER  MAINFRAME.  WHEN  THE  A/D  CARD 
IS  PROGRAMMED  TO  ZERO  VOLTS,  THE  OUTPUT  OF  THE  6224B  POWER 
SUPPLY  IS  -0.04S.  THEREFORE  AN  OFFSET  VOLTAGE  IS  ADDED  TO  THE 
PROGRAMMED  VOLTAGE.  THE  VOLTAGE  GAIN  IS  SET  TO  2.0. 

PARAMETERS - 


VH  -  HEATER  VOLTAGE  TO  EE  SET 

IMPORTANT  VARIABLES - 

GAIN  —  VOLTAGE  GAIN  =  6224B  VOLTAGE  /  PROGRAMMING  VOLTAGE 
OFFSET  —  NEGATIVE  CF  6224B  VOLTAGE  WHEN  ZERO  VOLTS  IS  PROGRAMMED 
RES  —  RESOLUTION  OF  A/D  CARD  =  5  MV  PER  BIT 
MAXVH  —  MAXIMUM  POSSIBLE  HEATER  VOLTAGE 

C 


REAL  VH,  GAIN,  OFFSET,  RES,  MAXVH 
INTEGER  BUF ( 2 ) ,  SLOT,  BITS 
C 

DATA  GAIN/2.0/,  OFFSET/0.0475/,  RES/C. 005/,  MAXVH/20.0/ 
DATA  BUF/170140B/,  SLOT/2O0OOB/ 

C 

IF  (VH.GT. MAXVH)  VH  =  MAXVH 
IF  (VH.LT.O.O)  VH  =  0.0 

BITS  =  IF IX  ((VH  +  OFFSET)  /  GAIN  /  EES  +  0.5) 

BUF ( 2 )  =  ICR  (SLOT,  BITS) 

C - EOT  ELEMENTS 

C  Bur ( 1 )  =  17014GB  —  CONTROL  WORD  DTE,  SYE  ON 

C  BUF ( 2 )  =  2XXXXP  —  A/D  CARD  DATA  WCRD 

C - MAKE  A  WRITE  WITH  HANDSHAKE  HECUEST 

C  ICGDS=2  CONTROL  WCRD=124B  BUFFER  L£NGTH=2  IrUNCsl 

CALL  EXEC  (2, 124B.EUF.2, 1 ) 

RETURN 

END 

END$ 
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C 

C*»*»*  RNS790802  -  60  SOURCE  FILE  :  RNTMQN 


•  • 

PROGRAM  TMON 

•  * 


THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  ACT  AS  AN  INTERFACE 
BETWEEN  THE  OPERATOR  AND  TEMPERATURE  CONTROL  PROGRAM. 

THE  MONITOR  COMMANDS  ARE  AS  FOLLOWS: 

T  —  CHANGE  TEMPERATURE  SET  POINT 
AND/OR  SAMPLING  PERIOD 
N  —  NO  CHANGE 

S  —  STOP  TEMPERATURE  CONTROL  PROGRAM 
AND  THIS  PROGRAM 

THE  SUBROUTINES  RDVM ,  TEMP  AND  MSTOP  ARE  INCLUDED  IN 
MULTI  PROGRAMMER  LIBRARY.  TO  LOAD  THEM  WITH  THIS  PROGRAM 
USE  THE  LOADER  COMMAND  SE.JMPLIB  TWO  TIMES. 

C 

INTEGER  PARAM( 5 ) ,  LU,  CLASS,  COMAND,  NAME ( 3 ) 

INTEGER  T,  N,  S,  DUMY1,  DUMY2 ,  1,  ERR,  CDELAY,  PDELAY 
REAL  MESSAG ( 2 ) ,  TCV,  TEMPER 
C 
C 

DATA  I/2HT  /,  N/2HN  /,  S/2HS  / 

C 

CALL  RMPAR ( PARAM ) 

LU  =  PARAMO) 

WRITE  (LU,  1 ) 

1  FORMAT  ("INPUT  THE  NAME  CF  THE  CONTROL  PROGRAM") 

READ  (LU , 2)  (NAME (I ) ,  1=1,3) 

2  FORMAT  (3A2) 

WRITE  (LU  ,  3  ) 

3  FORMAT  ("INPUT  REFERENCE  TEMPERATURE  AND  SAMPLING  PERIOD") 
READ  (LU , * )  ( MESSAG ( I ) ,  1=1,2) 

WRITE  ( L U  ,  4  ) 

4  FORMAT  ("AT  WHAT  INTERVAL  IN  SECONDS  DO  YOU  WANT  " 

$  "THE  TEMPERATURE  PRINTED  AT  THIS  TERMINAL  ?") 

READ  (LU,*)  PDELAY 
WRITE  (LU , 5 ) 

5  FORMAT  ("AT  WHAT  INTERVAL  IN  MINUTES  DO  YOU  WANT  " 

$  "TO  BE  PROMPTED  FOR  A  COMMAND  ?") 

PEAD  (LU,*)  CDELAY 
C 

C**##*  GET  A  CLASS  NUMBER  AND  PUT  MESSAG  IN  SYSTEM  E OFFER 
CLASS  =  0 

CALu  EXEC  ( 20,0, MESSAG, 4, DUMY 1 ,DUMY2, CLASS) 
r 
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C»***»  schedule  the  control  program  and  pass  it  the  class  nc. 

CALL  EXEC  (10,  NAME,  CLASS) 

GOTO  300 
C 

C*****  ENTER  THE  COMMAND  INPUT  LOOP 
100  WRITE  (LU , 6  ) 

6  FORMAT  ( "TMON  COMMAND  ?  ") 

READ  (LU ,2)  COMAND 

IF  (CCMAND  .EQ.  N)  GOTO  300 
IF  (COMAND  .NE.  T)  GOTO  500 
WRITE  (LU,3) 

READ  (LU, •)  (MESSAG(I),  1=1,2) 

WRITE  (LU , 4 ) 

READ  (LU,*)  PDELAY 
WRITE  (LU,5) 

READ  (LU,*)  CDELAY 

r 

:**•••  dq  A  CLASS  I/O  WRITE/READ 

CALL  EXEC  ( 20 , O.MESSAG , 4 , DUMY 1 , DUMY2 , CLASS) 

C 

C*»**»  READ  THE  TEMPERATURE  AND  PRINT  AT  THE  TERMINAL 

500  1=3 

UOO  CALL  RDVM(TCV) 

TEMPER  =  TEMP  (TCV) 

WRITE  (LU  ,  99  )  TEMPER 
99  FORMAT  ("TEMPERATURE  ="F10.3"  K" ) 

CALL  WAIT  (PDELAY,  2,  ERR) 

1  =  1  +  1 

IF  ( I*PDELAY/60  .LT.  CDELAY)  GOTO  4C0 
GOTO  100 
C 

500  IF  (CCMAND  .EQ.  S)  GOTO  600 
WRITS  (LU, 7 ) 

7  FORMAT  ( "ILLEGAL  COMMAND") 

GOTO  100 

C 

C***»*  KILL  ThE  CONTROL  PROGRAM  AND  RELEASE  THE  CLASS  NUMBER 
600  CALL  EXEC  (6,  NAME,  2) 

CALL  EXEC  (21,  CLASS,  MESSAG,  4) 

CALL  MSTOP 

END 

END$ 
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APPENDIX  B. 

Design  Algorithm  for  the  Discrete-Time  Controller 

The  design  procedures  covered  in  Section  2.1  may  be  put  in  the 
form  of  a  simple  design  algorithm  given  below: 

1.  Choose  the  damping  ratio  c  between  0  and  1.  A  damping  ratio 
near  1  gives  smaller  overshoots  of  the  step  response. 

2.  Find  ot  for  minimum  overshoot  with  the  given  value  of  c  from 
Equation  (2.11). 

3.  Choose  the  sampling  period  T. 

4.  Pick  a  reasonable  time  for  the  maximum  response,  t  .  This 

max 

time  should  be  greater  than  5J  of  the  plant  time  constant. 

5.  Find  the  angle  of  the  closed  loop  pole  p  from  Equation  (2.12). 

6.  Find  the  magnitude  of  the  closed  loop  pole  p1  from  Equation  (2.13). 

7.  Find  the  closed  loop  zero  z1  from  Equation  (2.14). 

8.  Find  the  closed  loop  constant  K  from  Equation  (2.4). 

9.  Find  the  controller  transfer  function  D(z)  from  Equation  (2.2). 

This  algorithm  is  ideal  for  execution  by  a  computer  design 

program.  A  program  which  performs  this  algorithm  is  listed  on  the 
following  pages.  A  sample  output  of  the  program  is  also  given  in 


this  appendix. 
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FTN4.L 

C 

RNS79063C  -  16  SOURCE  FILE  :  RNDSN 


C  *  * 

PROGRAM  DSN 

C  *  * 

C  »**t«»**»*#»#** 

r 

C  THIS  PROGRAM  PERFORMS  THE  DESIGN  ALGORITHM  FOR  A  DIGITAL 
C  CONTROLLER  WHICH  RESULTS  IN  A  CLOSED  LOOP  SYSTEM  WITH 
C  MINIMUM  STEP  RESPONSE  OVERSHOOT  FOR  A  GIVEN  VALUE  OF  THE 

C  DAMPING  RATIO,  ZETA.  THE  CLOSED  LOOP  SYSTEM  ALSO  POSSESES 

C  THE  MINIMUM  BANDWIDTH  NECESSARY  FOR  THE  PEAK  OF  THE  STEP 

C  RESPONSE  TO  BE  REACHED  AT  OR  BEFORE  A  SPECIFIED  TIME,  TM. 

C 

C  THE  CLOSED  LOOP  Z-TRANSFER  FUNCTION  IS  OF  THE  FORM: 

C 

C  K  *  (Z  -  Z1) 

C  M(Z)  =  - 

C  (Z  -  PI)  *  (Z  -  PI*) 

C 

C  WHERE  PI  AND  PI*  ARE  A  PAIR  OF  COMPLEX  CONJUGATE  POLES. 

C 

C  THE  Z-TRANSFER  FUNCTION  OF  THE  PLANT  IS  OF  THE  FORM: 

C 

C  GPK  *  (Z  -  GPZ 1 ) 

C  G?(Z)  - - - 

C  (Z  -  GPP  1 )  *  (Z  -  GPP2 ) 

c 

C  THE  RESULT  OF  THE  CALCULATION  IS  THE  TRANSFER  FUNCTION 
C  OF  THE  DIGITAL  CONTROLLER  : 

0 

0  NC( 1  )*Z**3  «■  NC ( 2 ) *Z**2  ♦  NC(3)*Z  +  NC(4) 

C  D(Z)  =  - - 

J  DC(  1  )*Z**3  *  DC(2 )*Z**2  t-  DC(3)*Z  ♦  DC(«) 


DESCRIPTION  OF  THE  VARIABLES 

0 '  -  CLOSED  LOOP  ZERO 

Pi  -  CLOSED  LOOP  COMPLEX  POLE 

<  .  CONSTANT  OF  THE  CLOSED  LOOP  TRANSFER  FUNCTION 

PHI  —  ARG ( P  1  ) 

WAGP 1  -  MAGNITUDE  OF  PI 
RE?i  --  »EAL  PART  OF  PI 
I  'PI  --  IMAGINARY  PART  CF  PI 

iPK  ---  INSTANT  OF  THE  PLANT  TRANSFER  FUNCTION 

,PZ1  --  ZERO  CF  THE  PLANT  TRANSFER  FUNCTION 

jPPI  --  Pr LE  ''y  THE  PLANT  TRANSFER  FUNCTION 

:?t>2  --  POLE  ir  THE  PLANT  TRANSFER  FUNCTION 

_E7A  —  DAMPING  RATIO 
CtTAC  -  OORT ( 1  -  ZETA**2) 


o  o 


r 


* 
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C  ALPHA  -  ARG(P1  -  Z1)  -  ARG(P1  -  1)  ♦  PI/2 

C  ALPHA  IS  CHOSEN  FOR  MINIMUM  OVERSHOOT 

C  T  -  SAMPLING  PERIOD  IN  SECONDS 

C  TM  — —  DESIRED  TIME  OF  PEAK  OF  STEP  RESPONSE 

C  NC  -  ARRAY  OF  NUMERATOR  COEFFICIENTS  OF  D(Z) 

C  DC  — —  ARRAY  OF  DENOMINATOR  COEFFICIENTS  OF  D(Z) 

C  A  -  TEMPORARY  VARIABLE 

C  9  -  TEMPORARY  VARIABLE 

C  BETA  —  TEMPORARY  VARIABLE 

C 

C 


DOUBLE  PRECISION  PI,  ZETA,  ZETA2,  T,  TM,  PHI,  BETA 
DOUBLE  PRECISION  MAGP1,  REP1 ,  IMP1,  Z 1 ,  K 
DOUBLE  PRECISION  GPK,  GPZ1,  GPP  1 ,  GPP2 
DOUBLE  PRECISION  NC(4),  DC(4),  A,  B 
C 

INTEGER  LU,  PARAM( 5 ) ,  I 
C 

DATA  PI/3. 141592654DO/ 

C 

CALL  RMPAR(PARAM) 

LU  =  PARAM( 1 ) 

C 

C  GET  THE  DESIGN  PARAMETERS  FROM  THE  USER 
WRITE  (LU, D 

1  FORMAT ( "  DIGITAL  CONTROLLER  DESIGN  PROGRAM"//) 

WRITE  (LU , 2) 

2  FORMAT ("  THE  Z-TRANSFER  FUNCTION  OF  THE  PLANT  IS  OF  THE  FORM  :"/) 
WRITE  (LU , 3 ) 

3  FORMAT (12X, "GPK  *  (Z  -  GPZD") 

WRITE  (LU , 4 ) 

4  FORMAT (  "  GPU)  - - ") 

WRITE  (LU , 5  ) 

5  FORMAT ( 9X ,  " ( Z  -  GPPD  *  (Z  -  GPP2) "/  ) 

WRITE  ( LU , 6 ) 

6  FORMATU  ENTER  THE  PARAMETERS:  GPK,  GPZ 1,  GPP1 ,  GPP 2") 

READ  (LU ,•)  GPK,  GPZ1,  GPP1,  GPP2 

100  WRITE  (LU,7 ) 

7  FORMATU"  ENTER  THE  DAMPING  RATIO,  ZETA,  BETWEEN  0  AND  1") 

READ  (LU,*)  ZETA 

IF  (ZETA.LE.O.ODO  .OR.  ZETA . GE. 1 . ODO )  GOTO  100 
200  WRITE  (LU , S ) 

8  FORMATU"  ENTER  THE  SAMPLING  PERIOD,  T,  IN  SECONDS") 

READ  (LU,*)  T 

IF  (T.LE.O.ODO)  GOTO  200 
300  WRITE  ( LU , 9 ) 

9  FORMATC/"  ENTER  THE  DESIRED  TIME  OF  THE  PEAK  STEP  RESPONSE,  TM" ) 
READ  (LU, *)  TM 

IF  (TM.LE.C.ODO)  GOTO  ?00 
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CALCULATE  ALPHA  FOR  MINIMUM  OVERSHOOT 
ZETA 2  =  DSGRT ( 1 . ODO  -  ZETA**2) 

ALPHA  =  DATN2( -ZETA ,  ZETA2 ) 


on  o  o  no  on 


c 

C  CHOOSE  PHI  FOR  MINIMUM  NECESSARY  BANDWIDTH 
PHI  =  PI  •  T  /  TM 

FIND  THE  CLOSED  LOOP  POLE,  PI 
MAGP1  =  DEXP ( -PHI  •  ZETA  /  ZETA2) 

REP1  =  MAG PI  *  DCOS(PHI) 

IMP1  =  MAG PI  *  DSIN(PHI) 

FIND  THE  CLOSED  LOOP  ZERO,  Z1 
BETA  =  ALPHA  +  DATN2HMP1,  REP1  -  1.0D0)  -  PI  /  2.0D0 
Z1  =  REP1  -  IMP1  •  DCOS(BETA)  /  DSIN(BETA) 

FIND  THE  CLOSED  LOOP  CONSTANT,  K 
K  =  ((I.ODO  -  REP1 ) #*2  ♦  IMP1 **2 )  /  (1.0D0  -  Z1) 

COMPUTE  THE  DENOMINATOR  COEFFICIENTS  OF  D(Z) 

A  =  2.0D)  *  REP1  ♦  K 
B  =  REP1 *#2  ♦  IMP  1 **2  ♦  K  •  Z1 
DC(U)  s  -B  *  GPZ1 
DC(3)  a  3  A  •  GPZ 1 
DC ( 2 )  =  -A  -  GPZ1 
DC( 1 )  =  1 . 0D0 

C 

C  COMPUTE  THE  NUMERATOR  COEFFICIENTS  OF  D(Z) 

A  s  GPP  1  +  GPP2 
B  :  GPP1  *  GPP2 
NC(1)  s  K  /  GPK 
NC(2 )  =  -NC(1)  •  (Z1  ♦  A) 

NC( 3 )  =  NC( 1  )  •  (Z1  •  A  ♦  B) 

NC(4)  =  -NC( 1 )  •  Z1  •  B 

c 

C  PRINT  THE  RESULTS 

WRITE  (LU.10) 

10  FORMAT( //"  THE  TRANSFER  FUNCTION  OF  TH:  DIGITAL  CONTROLLER  IS:"/) 
WRITE  (LU , 1 1 ) 

11  FORMAT (3X, "A0*Z**3  ♦  A1*Z**2  ♦  A2*Z  +  ,5") 

WRITE  (LU, 12) 

12  FOR  HAT  ("  D(Z)  - ") 

WRITE  (LU, 13) 

13  FOFMAT(3X,"EO*Z**3  +  B1«Z**2  +  B2*Z  B3"/) 

WRITE  (LU, 1U) 

14  FORMAT ( "  WHERE  :"/) 

DO  400  IsO, 3 

400  WRITE  (LU, 15)  I,  NC( 1+1) 

15  FOR MAT ( "  A" 1 1 "  =  "F10.6) 

DO  500  1=0,3 

500  WRITE  (IU,16)  I,  DC ( I ♦ 1 ) 

16  FOFMATC’  B"I1"  =  "F10.6) 

ENt 

ENT  $ 


A  sample  output  of  the  program  DSN  is  given  below: 


DIGITAL  CONTROLLER  DESIGN  PROGRAM 


THE  Z-TRANSFER  FUNCTION  OF  THE  PLANT  IS  OF  THE  FORM 

GPK  *  (Z  -  GPZ1 ) 

GPU)  . . 

(Z  -  GPP1 )  *  (Z  -  GPP2) 

ENTER  THE  PARAMETERS:  GPK,  GPZ1 ,  GPP1,  GPP2 
0.029964,  0.999600,  0.997603,  0.999660 

ENTER  THE  DAMPING  RATIO,  ZETA,  BETWEEN  0  AND  1 

0.8 

ENTER  THE  SAMPLING  PERIOD,  T,  IN  SECONDS 

2.0 

ENTER  THE  DESIRED  TIME  OF  THE  PEAK  STEP  RESPONSE,  TM 
90.0 


THE  TRANSFER  FUNCTION  OF  THE  DIGITAL  CONTROLLER  IS: 

A0*Z»«3  +  A 1  *Z**2  +  A2*Z  +  A3 

D(Z)  =  - . . 

B0*Z**3  +  S1»Z**2  +  B2*Z  +  B3 

WHERE  : 

AO  =  .212287 

A  1  =  -.224481 

A 2  =  -.186770 

A3  =  .198965 

BO  =  1.000000 

B 1  =  -2.823756 

E2  =  2.647582 

E  3  =  -.823826 
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elements.  The  principal  technical  mission  areas  are 
communication*,  electromagnetic  guidance  and  control,  sur- 
velllance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physic*  and  electronic  reliability,  maintainability  and 
compatibility. 
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